Maximum pull out force on DNA hybrids 

p. G. de Gennes 
College de France, 11 place Marcelin Berthelot 
75231 Paris Cedex 05, France 
E.mail: pgg@espci.fr 

February 2, 2008 



Abstract 

We discuss theoretically the force F between two colloidal particles, 
each of them carrying one single strand DNA. The two strands are com- 
plementary only on a finite sequence of {£) consecutive base pairs. We 
define an adjustment length x"^ (a few base pairs): in the adjustment 
regions near both ends of the paired region, the tension is still mainly on 
one single strand. But in the central part (for £ > >c~^) the two backbones 
are equally loaded. This leads to a rupture force Fc increasing linearly 
with £ for £>[<!, and saturating for £>c > 1. 



Abstract frangais 

Force maximum de tirage sur un ADN hybride 

Deux particules coUoidales, chacune portant une chaine greffee d'ADN, et les 
deux chaines pouvant s'hybrider, peuvent se ponter. On discute theoriquement 
la resistance en tension de ce pontage, lorsque I'hybrydation porte sur une 
sequence de £ bases consecutives. II apparait une longueur d'ajustement >f^^ 
(de I'ordre de quelques paires de base). Dans une region d'ajustement (de taille 
>c~^) pres de chaque extremite, la structure est distordue. Par contre, dans la 
region centrale (pour x£ >> 1) les deux chaines supportent des tensions egales, 
et la structure n'est pas perturbee. Ceci conduit a une force de rupture Fm qui 
augmente lineairement avec i pour >f£ < 1, et qui sature pour >ci > 1. 

Texte court en frangais 

Divers systemes de reconnaissance des acides nucleiques sont fondes sur des 
chaines a un brin, greffees sur une surface, et opposees a des chaines partielle- 
ment complementaires, greffees sur une autre surface -nanoparticule ou pointe 
d'un microscope de force [0, |Q, II existe une tension maximum Fc que 
pent supporter un pontage de ce genre. Dans la presente note, nous essayons 
de comprendre la forme de la relation Fc{£), et nous comparons les resultats a 
quelques donnees de la ref. [Q. 

Le modele utilise est un modele "d'echelle souple" decrit sur la figure |l|. 
Bien entendu, il faut imaginer cette echelle comme torsadee pour engendrer une 
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double helice, mais il semble que rhamiltonien de base (eq. |l|) garde un sens 
pour une forme plus realiste. 

Le resultat est une distortion, qui est importante pres des deux extremites de 
la zone appariee, et qui s'etend sur une longueur >c~^ pres des deux extremites. 
On ne gagne rien sur Fc a creer des structures hybrides de longueur £ » 
2>c^^. En comparant ces idees aux resultats de la reference [|l|, on arrive a une 
estimation ^ 6.2 paires de bases. 

1 Introduction 

Hybridisation between short DNA sequences, grafted on two surfaces, can lead 
to useful recognition systems. In particular, one can use an AFM tip §, ||, 
[^, measuring the force between tip and substrate for a single molecule. An 
obvious question then arises: what is the dependence of the rupture force Fc on 
the length of the hybridized portion? Is it useful to go to long lengths? 

In the present note, we try to provide a qualitative answer to this question, 
using a very crude "ladder model" for the hybridized portion. 

2 The ladder model 

We focus our attention on fig.|^, where the tension F is fed on one single strand 
at each end. A complete description of a distorted helix is, in principle, possible 
but complicated. We shall use here the simpler "ladder model" sketched on fig. 1. 
In the real world, the ladder is twisted into a double helix. But the distribution 
of forces should be qualitatively the same for both cases. 




Figure 1: Two particles linked by a single bridge 



Our starting point is a set of one dimensional displacement u„ and v„ for the 
two sides of the ladder: ie for the members of a base pair (n). m„ describes the 
{6e) portion (with n ranging from -£/2 to £/2) and v„ describes the conjugate 
portion {P"f). Un is different from v„, because the ladder is distorted by the 
force F. We postulate an elastic energy: 

oo i/2 f/2 

-Q {Un+l - Un) +2^^Q - V„) + 2^ {Un ~ V„) (1) 

-£/2 -oo -1/2 

The Q terms describe elongation of the backbone in one same piece: (0:7) or 
{5uj). The R terms come from the coupling between base pairs, and we expect 
R to be weaker than Q. 
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We shall supplement the elastic description of eq. |^ by a breaking condition: 
whenever the forces inside a base pair (n) are larger than a certain threshold 
/c, the bond will break. This corresponds to: 



i?|v„-u„|>/c (2) 
The equilibrium conditions derived from eq. (0) are: 



--^ — = Q (v„+i - 2v„ + v„_i) + R{un - v„) = (3) 

for all indices n in the interval {—i/2 ^ n ^ •^/2)- Outside of the interval, the R 
term drops out. We shall be concerned, in practice, with i values significantly 
larger than one, and go to the continuum limit: 



Q^ + i?(u-v) = (0) (4) 
There are similar equations for the other sequence: 



nil 

Q^+R{v-u) = (5) 
Adding (1) and (|), we find: 



Qi^(u + v)=0 (6) 
and this imposes a conservation of the total tension: 



Q-j^{u + y) = F ^ constant {n < e/2) 



' dn ' 

Un + v„ = nF/ Q 



(7) 

We then turn to a discussion of the difference Sn = u„— v„, which is ruled 



by: 

Q^-2RS^0 i\n\<e/2 (8) 

The solution is a combination of exponentials: for the problem at hand, the 
right combination is symmetric upon the exchange (n —^ —n): all the base pairs 
are distorted in the same direction: 

Sn — So cosh(>fn) (9) 
= 2R/Q (10) 



3 



If i? << Q, the "adjustment length" x ^ is larger than unity. 
The overall solution derived from eqs <M) and (lin) is: 



Un = nF/2Q + ^SoCoah{M:n) \ 
v„ — nF/2Q — i(5o cosh(><n) J 

The relation between Sq and F is derived from the boundary condition at n = 
£/2. Here, we must have: 



F = Q{Un - Un~l) + R{Un - V„) 

giving: 



F = Sa |0><sinh j + 2i?cosh (^>f- j | (12) 

The force on the last hydrogen bond {n ~ £/2) is i?(5„: when we reach the 
threshold /i, this corresponds to: 



MoCOsh(>^^/2) /i (13) 
Eq. (|l2|) then gives a global tension at threshold: 



Fc = 2/1 |>^-i tanh + 1 } (14) 

Two limits are of interest: 

a) short strands {k£ < 1), correspond to: 



F,^fi{e + 2) (15) 

Fc increases linearly with £. 

b) infinitely long strands: the rupture force reaches a maximum: 



Fc-^ F,n^2h{>i-^ + 1) (16) 

The force F^ is much larger than 2/i, because a number of base pairs 
work in parallel, near each end. 

Just below this maximum, we can write from eq. (n4) : 



Fc/2h = >c-\l-2e-'^') + l (17) 

All the base pairs in the adjustment regions (of length x^^) participate to 
the resistance, while the center portion is at rest. 

Note that, when the pairing breaks at both ends (n = ±^/2) the other ties 
inside {[n] < £/2), also break (since Fc is an increasing function of £). 
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3 Discussion 



1) Clearly, we do not gain strength by choosing a long pairing sequence: the 
optimal pairing length is 2k~^. This conclusion was achieved with a primitive 
"ladder model". But the existence and meaning of the length >r~^ are probably 
more general. 

2) The magnitude of >c~^ is unclear. If the leading coupling between strands 
is due to base pairing, this should be smaller than the covalent bonds, and x^^ 
should be larger than unity. But we can also think that the spring constant R 
in eq. (Q) is partly controled by the stacking of base pairs: bending the base 
plane by an angle ~ (0„— v„)/Z? (where D is the radius of the 2s DNA) may 
contribute another energy to i?, and reduce . 

If we start from the data of ref [Q, we find that the force Fc was roughly 
equal to 1.11 nanonewton for both I ~ 20 and i = 16: this corresponds to the 
saturation regime {Fc — Fm for Kt » 1). On the other hand, for t = 12, 
Fc = 0.83 nN. 

Using eqs (^, |l^), this leads us to >r^^ — 6.2 base pairs. Thus, at least for 
this particular example, the adjustment length seems to be much larger than 
unity. The corresponding force /i is 17 pN. 
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